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PCTEL-025 

METHOD AND APPARATUS FOR NONCOHERENT SIGNAL PROCESSING IN 

PILOTLESS WIRELESS SYSTEMS 

BACKGROUND OF THE INVENTION 

I. Field of the Invention 

The present invention relates to telecommunications. The invention more 
particularly relates to wireless telecommunications apparatus, systems and methods 
which implement data transmission via a plurality of telecommunication channels such as 
radio channels with variable parameters, including multipath wireless channels. More 
specifically, the invention relates to pilotless wireless systems with mobile transmitters 
and/or receivers, although it is not limited thereto. 



2. Stateof the Art 

Development of pilotless data transmission and signal processing is an important 
problem of wireless system design. Pilotless wireless systems provide the highest real 
data rate by utilizing system capacity exclusively for data transmission and providing all 
receiving functions without any accompanying pilot signals. 

One manner of implementing a pilotless system is using coherent demodulation 
based on a reference signal and extracting all necessary information from signals-bearing 
data. In this case, precise estimation and tracking of the carrier phase plays a critical role. 
A conventional approach to carrier phase estimation is described in John Proakis, "Digital 
Communications", McGraw Hill, Fourth Edition, 2000, Sections 6.2.4 -6.2.5 which is 
hereby incorporated herein in its entirety. A new method of carrier phase tracking was 
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1 proposed in co-owned U.S. Serial No. 10/628,943 filed July 29, 2003, and entitled 

2 "Pilotless, Wireless, Telecommunications Apparatus, Systems and Methods" which is 

3 also hereby incorporated by reference herein in its entirety. In that invention, phase 

4 tracking is based on reducing and averaging differential quadrature components of 

5 received symbols. 
6 

7 Coherent processing, based on carrier phase tracking, provides the maximum 

8 possible performance in channels with comparatively slow phase changing and 

9 comparatively high signal-to-noise ratio (SNR). If one of those conditions is not 

1 0 satisfied, the coherent demodulator loses its advantages. For example, when the bit error 

1 1 rate BER> 0.01, phase estimation becomes less precise. The SNR penalty depends on the 

12 method of phase tracking but it may reach 1-2 dB. However, much more performance 

1 3 loss can be caused by fast phase changing or, especially, phase jumping. Fast phase 

14 changing and phase jumping are typical phenomenon in communications involving 

1 5 mobile clients because in the mobile environment the multipath configuration may 

1 6 change instantly. This change can cause coherent detection degradation even at a 

1 7 comparatively high SNR because any phase tracking algorithm, based on symbol 

1 8 averaging, is not capable of instantly estimating the phase changing. As a result, the 

19 coherent receiver provides a long sequence (burst) of errors after phase jumping even at a 

20 comparatively high SNR. This is the reason why wireless system designers consider 

21 approaches other than coherent processing for pilotless systems. 
22 

23 For multipath radio channels in a mobile environment, a promising manner of 

24 implementing a pilotless system involves the utilization of noncoherent signal processing 

25 which does not need any information about the initial phase of the reference signal and 

26 consequently does not require phase tracking at all. In the case of phase modulation, 
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1 noncoherent processing can be only used in combination with Differential Phase Shift 

2 Keying (DPSK) in contrast to coherent processing which can be used with both DESK 

3 and PSK. As is emphasized in Dariush Divsalar, and Marvin Simon, "Multiple-Symbol 

4 Differential Detection of MPSK", IEEE Transactions on Conmiunications, vol.38, N3, 

5 March 1990 (which is hereby incorporated by reference in its entirety), for multipath 

6 channels with fast phase changing, a DPSK-with-noncoherent-processing scheme is the 

7 only way to provide robust data transmission. 
8 

9 Conventional noncoherent DPSK processing includes a two-symbol interval for 

10 making a decision. With a two-symbol interval, if a phase jump occurs during the 

1 1 symbol interval, only one error will take place. However, two-symbol non-coherent 

12 DPSK suffers from a performance penalty when compared to ideal coherent DPSK. For 

1 3 example, at BER=10"^ the penalty is about 0.75 dB for DBPSK (Differential Binary 

14 PSK), and about 2.2 dB for DQPSK (Differential Quadrature PSK). For M-ary DPSK 

1 5 (DMPSK) the penalty increases with increasing M. 
16 

17 A well-known way to mitigate the performance loss of the non-coherent receiver 

18 and preserving its advantages is to increase the interval of the non-coherent processing. 

1 9 This approach, known as multisymbol non-coherent detection (or processing) was 

20 considered in the previously incorporated article by Dariush Divsalar, and Marvin Simon. 

21 The conventional approach to multisymbol processing is based on allowing the 

22 observation interval over which symbol decisions are made to be more than a two-symbol 

23 interval while at the same time making a joint decision on several symbols 

24 simultaneously as opposed to symbol-by-symbol decisions. An obvious disadvantage of 

25 the joint decision on several symbols is the additional symbol delay which results in the 
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1 receiver. In addition, the joint decision procedure is undesirable for some decoding 

2 algorithms. 
3 

4 SUMMARY OF THE INVENTION 

5 

6 It is therefore an object of the invention to provide pilotless teleconmiunication 

7 systems and associated methods and apparatus, which transmit data by means of DPSK 

8 modulation of carriers and which provide all receiving functions based on multisymbol 

9 noncoherent processing of carriers. 
10 



11 It is another object of the invention to provide general methods and apparatus for 

12 multisymbol noncoherent processing in single-carrier and multi-carrier systems with 

13 symbol-by-symbol decisions. 
14 

15 It is a further object of the invention to provide relatively simply implementable 

1 6 algorithms for multisymbol noncoherent processing of M-ary DPSK modulated signals. 
17 

18 It is an additional object of the invention to provide simple iterative algorithms of 

19 multisymbol noncoherent processing of DPSK modulated signals, when each next 

20 symbol may be simply combined with a result of previous symbol processing(s). 
21 

22 Another object of the invention is to provide pilotless wireless telecommunication 

23 systems, apparatus, and methods which provide all receiving functions without pilot 

24 signals by means of multisymbol noncoherent processing with symbol-by-symbol 

25 decisions which are based on iterative algorithms. 



26 
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1 In accord with the above objects which will be discussed in more detail below, the 

2 present invention provides methods, and apparatus for the realization of pilotless 

3 teleconmiunication systems, which use the entire system capacity exclusively for data 

4 transmission, and which provide all receiving functions based on an iterative multisymbol 

5 noncoherent processing of carriers with symbol-by-symbol decision-making. 
6 

7 According to the invention, a method or algorithm is provided for 3-symbol 

8 noncoherent processing of M-ary DPSK (DMPSK) modulation with an arbitrary set of 

9 phase differences. The preferred algorithm is a three-step procedure. The first two steps 

10 are iterative, with the first step being the calculation of two-symbol quadrature 

1 1 components, and the second step being the calculation of three-symbol quadrature 

12 components of final three-symbol vectors. The lengths of the vectors are then compared 

13 in the third step to find the maximum. The algorithm is preferably implemented with a 3- 

14 symbol DPSK-signal noncoherent demodulator which utilizes intersymbol processors, 

1 5 memory registers and a decision block. 
16 

17 The preferred algorithm of the invention has a desirable iterative structure in that 

18 the 3-symbol quadrature components are sums of the previous 2-symbol components and 

1 9 transformed quadrature components of the current symbol. The preferred algorithm of 

20 the invention provides the ability to make a symbol-by-symbol decision and/or a joint 

21 decision on two symbols simultaneously. 
22 

23 According to another aspect of the invention, a simplified algorithm for 3-symbol 

24 noncoherent processing of DQPSK signals is provided. 
25 
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1 According to a further aspect of the invention, a general algorithm for N-symbol 

2 noncoherent processing of M-ary DPSK signals is provided. The general algorithm is 

3 carried out with (N-1) recurrent steps (iterations) plus a decision step. Each iterative step 

4 includes a simple trigonometrical transformation of quadrature components of the current 

5 symbol and a sum of the transforms with the results of the previous step. A final N- 

6 symbol based decision regarding the current transmitted symbol corresponds to the vector 

7 of maximum length, calculated after the (N-l)-th step of the recurrent procedure. The 

8 general algorithm is optionally implemented with (N-1) intersymbol processors, (N-2) 

9 pairs of memory blocks for saving results of the intersymbol processors, a decision block, 

10 and two shift registers for quadrature components of the received signal. Alternatively, a 

1 1 single processor may be used for all functions of the intersymbol processors, and one or 

12 more memory blocks may used to implement the memory blocks and shift registers. 
13 

14 Additional objects and advantages of the invention will become apparent to those 

15 skilled in the art upon reference to the detailed description taken in conjunction with the 

1 6 provided figures. 
17 
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1 BRIEF DESCRIPTION OF THE DRAWINGS 
2 

3 Figure 1 is a high level block diagram of a single carrier receiver incorporating a 

4 multisymbol noncoherent demodulator according to the invention. 
5 

6 Figure 2 is a high level block diagram of a multicarrier receiver incorporating a 

7 multisymbol noncoherent demodulator according to the invention. 
8 

9 Figure 3 is a schematic diagram of a 3-symbol noncoherent demodulator of DPSK 

1 0 signals according to the invention. 
11 

12 Figure 4 is a schematic diagram of the intersymbol processor and memory 

1 3 registers of the 3-symbol noncoherent demodulator of Fig. 3. 
14 

15 Figure 5 is a schematic diagram of the decision block of the 3-symbol 

1 6 noncoherent demodulator of Fig. 3 . 
17 

18 Figure 6 is a schematic diagram of an N-symbol noncoherent demodulator of 

1 9 DPSK signals according to a second embodiment of the invention. 
20 

21 Figure 7 is a high level schematic diagram of an alternative second embodiment 

22 of the invention which utilizes a high speed intersymbol processor in an N-symbol 

23 noncoherent demodulator of DPSK signals. 
24 

25 Figure 8 is a graph showing results of a stochastic simulation of the algorithms of 

26 the invention, with the bit error rate (BER) plotted as a function of signal-to-noise ratio 
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1 (SNR) for DBPSK modulation and noncoherent processing with 2, 3, 4, and 5-symbol 

2 intervals. 
3 

4 Figure 9 is a graph showing results of a stochastic simulation of the algorithms of 

5 the invention, with the bit error rate (BER) plotted as a function of signal-to-noise ratio 

6 (SNR) for DQPSK modulation and noncoherent processing with 2, 3, 4, and 5-symbol 

7 intervals, 
8 

9 

1 0 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

11 

12 Turning now to Fig. 1, a block diagram of a single carrier receiver 10 is shown. 

13 Receiver 10 includes a channel interface 12, an in-phase/quadrature (I/Q) detector 14, a 

14 multisymbol noncoherent demodulator of the invention 20, and a decoder 22. In the case 

15 of a wireless system, the channel interface 12 typically includes an antenna, a low noise 

1 6 amplifier, and possibly a signal transformer which transforms the signal to an intermediate 

17 frequency. The I/Q detector 14 typically includes two multipliers 23a, 23b, two low pass 

1 8 filters 24a, 24b, a local oscillator 25, and a phase rotator 26. As is well known in the art, 

1 9 the I/Q detector takes the signal (possibly, the intermediate frequency signal) from the 

20 channel interface and provides quadrature components X, Y to the demodulator 20. 

21 Determinations of the demodulator 20 are provided to the decoder 22. According to the 

22 invention, different implementations of channel interfaces, I/Q detectors and decoders 

23 may be utilized as are known in the art. 
24 

25 A block diagram of an OFDM (orthogonal frequency division multiplexed) 

26 multicarrier receiver 10' is shown in Fig. 2. Receiver 10' includes a channel interface 12, 
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1 an I/Q detector 14, a fast Fourier transformer (FFT) 15, a multisymbol noncoherent 

2 demodulator of the invention 20', and a decoder 22\ The channel interface 12 preferably 

3 includes the same elements as in the single carrier system of Fig. 1. The FFT 15 provides 

4 quadrature components of all carriers for the received signal. The output of the FFT is a 

5 set of quadrature components X^n and Y^n for carrier numbers (cn)=l,2,...L where L is the 

6 number of carriers of the multicarrier system. The multisymbol noncoherent demodulator 

7 20' processes each pair of quadrature components in parallel or in a consecutive manner 

8 depending on desired implementation. The output of the demodulator 20' is provided to 

9 decoder 22*. Again it will be appreciated that different implementations of channel 

1 0 interfaces, I/Q detectors, Fourier transformers (or equivalents), and decoders may be used 

11 as are known in the art. 
12 

1 3 Prior to turning to the figures showing the noncoherent demodulators of the 

14 invention, a discussion of the theoretical underpinnings of the demodulators of the 

15 invention is useful. First, it is desirable to consider 2-symbol and 3-symbol noncoherent 

1 6 processing of M-ary DPSK (DMPSK) modulation with a set of phase differences 

17 A(k); k=l,2,...,M. (1) 
18 

19 In particular, define X(q) and Y(q) as the quadrature components (real and imaginary 

20 parts) of the received signal at the q-th symbol interval, where q=l , 2, 3, . . . It should be 

21 noted that quadrature components X(q) and Y(q) are results of the convolution of the 

22 received signal and the reference signal, or results of a Discrete Fourier Transform (DFT) 

23 of the received signal based on a certain reference signal. For purposes of this discussion 

24 it is assumed that the frequency of the reference signal is equal or "close enough" to the 

25 frequency of the received signal. This means that the system preferably provides some 

26 sort of frequency offset compensation, which is one of the conventional receiver 
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1 functions. Frequency olfset compensation allows a receiver to mitigate SNR degradation 

2 in the case of single carrier systems and avoid intercarrier interferences in the case of 

3 multicarrier OFDM systems. 
4 

5 With these circumstances, the decision-making procedure for conventional 2- 

6 symbol processing comprises two operations. The first operation is the calculation of the 

7 2-symbol quadrature components VXj and VY2 for all versions of the transmitted phase 

8 differences A(k) according to: 
9 

10 VX2 (kj) = X(q-l) + X(q)cos[A(kJ] + Y(q)sin[A(k,)], (2a) 

1 1 YY, (k,) = Y(q-l) + Y(q)cos[A(k,)] - X(q)sin[A(k,)], (2b) 

12 whereq=2, 3, and ki=l, 2, M. 
13 

14 In equations (2a) and (2b), the consequent (current) symbol is reduced to the most likely 

15 preceding symbols. It is possible to use the inverse transformation: reduction of the 

16 preceding symbol into the most likely consequent symbols. Which option is utilized 

17 depends on the direction of the phase difference reading, but the final result of the 

18 decision making will be the same. In the following, reduction of the consequent symbol 

1 9 into the possible preceding symbols is utilized. 
20 

21 The second operation for conventional 2-symbol processing is finding the index 

22 "k," corresponding to the maximum vector length: 
23 

24 max{ [VX^ (k,)]' + [VY^ (k^)]^}; k,=l, 2, M. (3) 

25 
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1 It is noted that M sets of quadrature values are calculated in the first operation, and the M 

2 values in the square brackets of (3) are compared to find the maximum. 
3 

4 According to the invention, a symbol-by-symbol decision-making algorithm for a 

5 3-symbol processing interval starts with the determination of the quadrature components 

6 VX2 and VYj. These components in the 3-symbol case, for the q-th, (q- l)-th and (q-2)-th 

7 symbols are equal to 
8 

9 VX2 (k,) = X(q-2) + X(q-l)cos[A(k,)] + Y(q-l)sin[A(k,)], (4a) 

10 VYj (k,) = Y(q-2) + Y(q-l)cos[A(k,)] - X(q-l)sin[A(k,)]. (4b) 
11 

1 2 Variables UX3 and UY3 can now be defined corresponding to the transformation of the 

1 3 quadrature components X(q) and Y(q) of the current q-th symbol for each possible 

14 combination of phase differences A(k,) between the (q-2)-th and (q-l)-th symbols, and 

1 5 A(k2) between the (q- 1 )-th and q-th symbols according to: 

16 
17 

1 8 UX3(k„k2) = X(q)cos[A(k,)+ A(k2)] + Y(q)sin[A(k,) + ^(k^)] (5a) 

19 UY3(k„k2)= Y(q)cos[A(k,) + AOc^)] - X(q)sin[A(k,) + m,)]. (5b) 
20 

21 Now, to find the 3-symbol quadrature components VXj and VY3 , transforms .(5) are 

22 added to the previously calculated values (4): 
23 

24 VX3(k„ kj) = VX2(k,) + UX3(k„ kj); (6a) 

25 VY3(k„ k^) = VY2(k,) + UY3(k„ k^). (6b) 
26 
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1 Finally, a decision is made by finding a maximum: 
2 

3 Max { [VX3(k., k^)]^ + [VY.iK M'}* (6c) 

4 wherek,=l,2,...,M, and k2=l,2,...,M. 
5 

6 So, the algorithm for 3-symbol processing interval can be described as a three step 

7 procedure: the first step is calculating 2-symbol quadrature components VX2and VYj 

8 according to equations (4a) and (4b), and the second step is calculating 3-symbol 

9 quadrature components VX3 and V Y3 according to equations (6a) and (6b). The third 
10 step is making a decision by finding the maximum value according to equation (6c). 
11 

12 According to equations (6a) and (6b), M components of VX3 and M components 

13 of VY3 are calculated for each component VX2 and VYj, i.e., of pairs (VX3, VY3) can 

14 be determined. According to the invention, the pairs of vector lengths are then 

1 5 compared to find the maximum according to equation (6c). For example, in the case of 

16 DQPSK (as discussed below), sixteen pairs of (VX3, VY3) are calculated and compared. 
17 

18 Two additional notes should be made regarding equations (6a) - (6c). First, 

19 equations (6a) and (6b) have a desirable iterative structure in that the 3-symbol 

20 quadrature components VX3, VY3 are the sums of the previous 2-symbol components 

21 VX2, VY2 and the transformed quadrature components of the current symbol. Second, 

22 the decision algorithm (6c) may be used for both a symbol-by-symbol decision and/or a 

23 joint decision on two symbols simultaneously. In the case of a symbol-by-symbol 

24 decision, maximization in (6c) is provided through indexes k2, corresponding to phase 

25 differences between the (q-l)-th and q-th symbols. In the case of the joint symbol 
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1 decision, maximization in (6c) is provided through the whole set of indexes k, and kj, 

2 corresponding to two sets of phase differences: i.e., the phase difference between the 

3 (q-2)-th and the (q- l)-th symbols, and the phase difference between the (q-l)-th and the 

4 q-th symbols. 
5 



6 To reiterate, according to the invention, an algorithm for a 3-symbol noncoherent 

7 processing of DMPSK telecommunications signals can be described as follows: 

8 a) Three consecutive symbols with indexes q, (q-1) and (q-2) are used to participate 

9 in making a decision relative to the current (q-th) symbol, and the corresponding 

10 quadrature components of those signals X(q) and Y(q), X(q-l) and Y(q-l), X(q-2) 

1 1 and Y(q-2) are saved; 

12 b) Quadrature components X(q- 1 ), Y(q- 1 ) and X(q-2), Y(q-2) are transformed into a 

1 3 set of 2-symbol quadrature components VXj and VYj according to equations (4a) 

1 4 and (4b) using each possible phase difference A(ki), where ki= 1 , 2, . . . , M; 

15 c) Quadrature components X(q) and Y(q) of the current symbol are transformed into 

16 a set of transforms UX3(ki, kj) and UY3(ki, kj) according to equations (5a) and 

1 7 (5b) for each possible combination of phase differences A(k,) and A(k2), where 

18 k,=l,2,...,M and k2=l, 2,..., M; 

19 d) Sets of 3-symbol quadrature components VX3(ki, k2) and VY3(k|, kz) are 

20 combined by adding the results of b) and c) according to equations (6a) and (6b); 

21 e) A final decision relative to the current (q-th) symbol is the phase difference of the 

22 set of differences A(k2), where k2=l, 2,. . ., M, corresponding to the maximum 

23 length of vectors [VX3(ki, k2),VY3(k,, k2)], according to equation (6c). 
24 

25 Fig. 3 shows a schematic diagram of the 3-symbol noncoherent demodulator of 

26 DPSK signals 20a (corresponding to demodulator 20 of Fig. 1, or 20' of Fig. 2) according 
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1 to the invention, operating according to a) through e) above. In particular, the 

2 demodulator 20a of the invention includes four identical delay elements D„ Dj, Dj and D4 

3 with delay times equal to symbol length T, two memories Rj and Rj, two intersymboi 

4 processors IP, and IP2, and 3-symbol decision block Decision3. 
5 

6 The demodulator 20a operates as follows. Quadrature components X(q) and Y(q) 

7 of the q-th received symbol are directly fed to intersymboi processor IP2 and through 

8 delay elements D1-D4 to intersymboi processors IP,. In IP,, two adjacent symbols 

9 (X(q-l),Y(q-l)) and (X(q-2),Y(q-2)), preceding the current symbol (and found in the 

10 delay elements Di, Dj, D3 and D4), are processed according to equations (4a) and (4b) or 

1 1 according to Table 1 (below) for DQPSK (as described hereinafter) to generate M pairs 

12 of numbers VXj and VYj which are stored in memories R, and Rj, respectively. In the 

1 3 case of DQPSK (as described below), the memories R, and Rj store four numbers for 

1 4 VX2 and four numbers for V Yj. The numbers VXj and V Yj from memories R, and Rj 

15 together with current symbol (X(q),Y(q)) are fed to IP2, where they are processed 

16 according to equations (5a), (5b) and (6a), and (6b) (or according to Tables 5a and 5b for 

17 DQPSK as described below). Intersymboi processor IP2 generates numbers VX3 and VY3 

18 which are fed to the decision block Decisions, where a decision is made according to 

19 equation (6c). In the case of DQPSK, as described below, sixteen pairs of numbers 

20 (VX3, VY3) participate in the decision making procedure. 
21 

22 According to one aspect of the invention, equations (6a) and (6b) may be 

23 considerably simplified in certain circumstances. In particular, in DQPSK modulation 

24 systems, a typical set of phase differences are: 
25 

26 A(k)= (0, Ji/2, Jt, 3jr/2). (7) 
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1 With that limited set of phase differences, values VXj and VYj calculated according to 

2 equations (4a) and (4b) can be reduced as suggested in Table 1. 
3 

4 Table 1: Values VXj and VYj according to equations (4a) and (4b) 



k, 


1 


2 


3 


4 


A(k,) 


0 


31/2 


n 






X(q-2)+X(q-l) 


X(q-2)+Y(q-l) 


X(q-2)-X(q-l) 


X(q-2)-Y(q-l) 


VY, 


Y(q-2)+Y(q-l) 


Y(q-2)- X(q-l) 


Y(q-2)-Y(q-l) 


Y(q-2)+X(q-l) 



5 
6 

7 From Table 1 it can be seen that the calculation of values VXj and VYj are reduced to 

8 simple combinations (summation or subtraction) of quadrature components of the 

9 adjacent received symbols. 
10 

1 1 Table 2 contains possible values of the sum (A(k,)+ A(k2))modulo 2k found in the 

12 square brackets of equations (5a) and (5b) for the DQPSK demodulator with phase 

13 differences (7). 
14 

1 5 Table 2: All combinations of [A(k,)+ A(k2)] modulo 2jt for DQPSK. 



\k2 
k^ 


1 


2 


3 


4 


1 


0 


Td2 


n 


3jt/2 


2 


n/2 




3jt/2 


0 


3 


n 


3jt/2 


0 


nl2 


4 


3JI/2 


0 


Jt/2 





16 
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1 
2 
3 
4 
5 



6 
7 
8 



9 
10 
11 
12 
13 
14 



Continuing with equations (5a) and (5b) for DQPSK, Tables 3a and 3b contain 
sine and cosine functions of the combinations indicated in Table 2: 

Table 3a: Sin[A(k,)+ A(k2)] 



k, 


1 


2 


3 


4 


1 


0 


1 


0 


-1 


2 


1 


0 


-1 


0 


3 


0 


-1 


0 


1 


4 


-1 


0 


1 


0 


3b: Cos[A(k,)+ A(k2)] 




1 


2 


3 


4 


1 


1 


0 


-1 


0 


2 


0 


-1 


0 


1 


3 


-1 


0 


1 


0 


4 


0 


1 


0 


-1 



With the sine and cosine functions determined. Tables 4a and 4b contain values 
for the transforms UXjCk,, kj) and UY3(ki, kj) of equations (5a) and (5b). 
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Table 4a: UXjOc,, kj) = X(q)cos[A(k,) + ACk^)] + Y(q)sin[A(k,) + A(kj)] 



K 


1 


2 


3 


4 


1 


X(q) 


Y(q) 


-X(q) 


-Y(q) 


2 


Y(q) 


-X(q) 


-Y(q) 


X(q) 


3 


-X(q) 


-Y(q) 


X(q) 


Y(q) 


4 


-Y(q) 


X(q) 


Y(q) 


-X(q) 



Table 4b: UYjCk,, k^) = Y(q)cos[A(k,) + ACk^)] - X(q)sin[A(k,) + AOcj)] 





1 


2 


3 


4 


k,\ 










1 


Y(q) 


-X(q) 


-Y(q) 


X(q) 


2 


-X(q) 


-Y(q) 


X(q) 


Y(q) 


3 


-Y(q) 


X(q) 


Y(q) 


-X(q) 


4 


X(q) 


Y(q) 


-X(q) 


-Y(q) 



It should be noted that each row in Tables 2-4 is a cyclic shift of the previous row. This 
fact gives additional opportunity for simplification of the computation procedure as it is 
necessary to store only four numbers for the computation of the final sixteen numbers. 

Finally, using the results from Tables 4a and 4b, Tables 5a and 5b contain the 3- 
symbol quadrature components VX3 and VY3, calculated according to equations (6a) and 
(6b). 



PCTEL-025 



1 Table 5a: VXj (k„ k,) = VXj (k,)+ UXjOc,, kj) 





1 


2 


3 


4 


1 


VX2(l)+X(q) 


VX2(l)+Y(q) 


VX2(l)-X(q) 


VX,(1) -Y(q) 


2 


VXj(2)+Y(q) 


VX2(2)-X(q) 


VX,(2) -Y(q) 


VXj(2)+X(q) 


3 


VXjO) -X(q) 


VX2(3)-Y(q) 


VX,(3)+X(q) 


VX2(3)+Y(q) 


4 


VX2(4) -Y(q) 


VXj(4)+X(q) 


VX,(4)+Y(q) 


VXj(4) -X(q) 



2 
3 

4 Table 5b: VY3(k„ k^) = VYjCk,) + UYjCk,, k^) 





1 


2 


3 


4 


1 


VY2(1)+Y(q) 


VY^d) -X(q) 


VY^d) -Y(q) 


VYa(l)+X(q) 


2 


VY2(2)-X(q) 


VY2(2) -Y(q) 


VY2(2)+X(q) 


VY2(2)+Y(q) 


3 


VY^O) -Y(q) 


VY2(2) +X(q) 


VY2(2)+Y(q) 


VY2(2) -X(q) 


4 


VY2(4)+X(q) 


VY2(4)+Y(q) 


VY,(4) -X(q) 


VY2(4) -Y(q) 



5 

6 The components of Tables 5a and 5b provide the basis for making a final decision 

7 according to algorithm (6c). In the case of DQPSK there are four possible phase 

8 differences, (i.e., the differences set forth in (7)), and for purpose of simplification we can 

9 designate the corresponding decisions with numbers 1, 2, 3 and 4. These numbers for the 

1 0 third symbol of 3-symbol DQPSK processing are shown in Table 6. The decision 

1 1 corresponds to the location of the maximum according to equation (6c). This means that 

12 if, for example, the maximum of [VX3(ki, kj)]^ + [VY3(k„ k^)? is achieved at k,=2 and 

13 k2=3, then the decision is 3. So, for a one-symbol decision it is sufficient to determine a 

14 column number in Table 6, in which the maximum of equation (6c) takes place. 
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1 Table 6. Decisions for the third symbol of 3-symboI processing of DQPSK. 





1 


2 


3 


4 












1 


1 


2 


3 


4 


2 


1 


2 


3 


4 


3 


1 


2 


3 


4 


4 


1 


2 


3 


4 



2 

3 Joint decisions for the second and third symbols of the 3-symbol DQPSK 

4 processing are shown in Table 7, As in the previous case, a decision corresponds to 

5 location of the maximum according to equation (6c). Thus, for example, if the maximum 

6 of [VXjCki, kj)]^ + [VY3(k„ kj)]^ is achieved at ki=3 and k2=4, then the decision is 3 for 

7 the second symbol, and the decision is 4 for the third symbol. Generally, all sixteen 

8 components [YXjCkj, kj)]^ + [VYjCki, kj)]^ should be compared to find the maximum. 
9 

10 Table 7. Joint decisions for the second and third symbols of 3-symbol 

1 1 processing of DQPSK. 
12 





1 


2 


3 


4 


1 


1.1 


1,2 


1,3 


1,4 


2 


2,1 


2,2 


2,3 


2,4 


3 


3.1 


3,2 


3,3 


3,4 


4 


4,1 


4,2 


4,3 


4,4 



13 
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1 Turning back to the schematic diagram in Fig. 3, it should be appreciated that in 

2 the case of DQPSK, the intersymboi processor IP, provides the computations of Table 1, 

3 and the intersymboi processor IP2 provides the computations of Tables 5. 
4 

5 Fig. 4 shows additional details of the intersymboi processor IP2 and memories R, 

6 and for the 3-symbol noncoherent demodulator of DQPSK signals. The intersymboi 

7 processor IP2, shown within dotted lines, includes a four-cell circular shift register (CSR) 

8 and two adders Ai and A2. Memories Ri and R2are four-cell shift registers. As seen in 

9 Fig. 4, there are four signals received at the input to the memories R, and Rj and to 

10 intersymboi processor IP2 including: the current quadrature components X(q) and Y(q) of 

11 the received signal, and numbers VXj and V Y2 from the output of the previous 

1 2 intersymboi processor IPp Numbers X(q) and Y(q) are saved in the four cells of the CSF 

13 as positive and negative values as indicated in Fig. 4 and as required by Tables 4a and 4b. 

14 In addition, four numbers VXj are saved in register R| and four numbers of VY2 are 

1 5 saved in register R2, all of which are output from intersymboi processor IPj. 
16 

17 The calculation procedure in the intersymboi processor IP2 consists of four 

1 8 cycles, each cycle having four subcycles. During the four subcycles of a first cycle, the 

19 first components of VX2 and VY2 are sequentially summed in the adders Ai and A2 with 

20 all four elements of the CSR (with the numbers in the CSR circularly shifting or rotating 

21 between each subcycle). The first cycle generates the first rows of Tables 5a and 5b. As 

22 seen in Fig. 4, the addition of VX2 and VY2 with numbers in the CSR is done in parallel, 

23 with VX2 being added to the number in the last cell of the CSR, and V Y2 being added to 

24 the next-to-last cell of the CSR. After the first cycle, the numbers in the shift registers R, 

25 and R2 are shifted such that the next component of VX2 and V Yj may be sequentially 

26 sunmied in the adders A, and Aj with all four elements of the CSR to provide the next 
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1 line of Tables 5a and 5b. The process continues until all four components of VXj and 

2 VY2 stored in registers Ri and R2 are added with all four elements of the CSR according 

3 to Tables 5a and 5b. After each cycle, numbers in shift registers R| and Rj are shifted, 

4 setting free space for the next symbol data. 
5 

6 The results of the interprocessor IPj is to provide sixteen sets of values (VXj.VYj) 

7 for the decision block. Fig. 5 shows a schematic diagram of the decision block. Two 

8 numbers at a time from the output of the intersymbol processor IP2 are fed to squarers or 

9 multipliers M| and Mj, the outputs of which are then combined in an adder A3 in order to 

10 implement relationship (6c). The resulting sequence of numbers from the output of 

1 1 the adder A3 is fed to a logic element 30 which compares each incoming number to a 

12 variable V which is initialized as V=0, If at comparator 32, the current number > V, 

1 3 then V is replaced by Wj in a memory 34 of variable V. In total, sixteen numbers Wj are 

14 compared in the case of the 3-symbol processing of DQPSK signals to find the maximum 

1 5 value. If a decision for only the third symbol is desired, the index "i" of number Wj at the 

1 6 output of the comparator 32 can be transformed at 36 into j=i mod4 such that each new 

1 7 number "j" replaces the previous one in a memory 38 of variable j. According to Table 6, 

1 8 the final value of variable "j" is the desired decision for the third symbol. However, for a 

19 joint two-symbol decision (i.e., decisions for the second and third symbols), logic unit 30 

20 should operate according to Table 7, with index "i" defining both k| and in order to 

21 provide a decision for a single element of Table 7 with maximum value . Thus, for 

22 example, if W7 proved to be a maximum, k,=2 and k2=3, the decision would be 2 for the 

23 second symbol and 3 for the third symbol Similarly, if W9 proved to be a maximum, 

24 kj=3 and k2=l, the decision would be 3 for the second symbol and 1 for the third symbol. 
25 
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1 Given the above, it will be appreciated by those slcilled in the art that the 

2 invention can be generalized to provide a symbol-by*symbol decision making algorithm 

3 for M-ary DPSK at an N-symbol processing interval. In particular, consider DPSK 

4 modulation with an arbitrary set of phase differences. The algorithm of N-symbol 

5 processing of the q-th received symbol X(q),Y(q) is carried out with (N-1) recurrent steps 

6 (iterations). 
7 

8 Let VXi^.|(k„ kj, . . M kj) and VYi+i(k„ kj, . . be quadrature components at the 

9 i-th step, and UXi+i(k,, kz, . . kj) and UYi+i(ki, kj, . . ., kj) be transforms of the received 

10 symbol at the i-th step, where i=l, 2,.. ., N-1. These components and transforms are 

1 1 functions of integer variables k,, kj, . . Iq, running from 1 to M. Then, transforms 

1 2 UXi+,(ki, kj, . . . , kj) and UYi+,(k„ kj, . . . , k^) of the received symbol at the i-th step can be 

13 calculated according to the following which is similar to equations (5a) and (5b): 
14 

15 UXi^,(ki, kj, ki) = X(q-N+i+l)cos[A(k,)+...+ A(ki)] + Y(q-N+i+l)sin[A(k,)+...+ A(ki)] (8a) 

16 UYi^i(k„ k2. ki) = Y(q-N+i+l)cos[A(k,)+...+ A(ki)] - X(q-N+i+l)sin[A(k,)+...+ A(ki)], (8b) 

1 7 where q is a symbol number which is currently under the decision making procedure. 
18 

19 Now the i-th step (iteration) of the desired recurrent algorithm can be represented 

20 with the following simple formula: 

21 

22 VXi,i(k„k2, ...,ki)= VXi(ki,k2, ...,k,,) + UXi,,(k„k2, ...,ki), (9a) 

23 VY;,,(k„ k^, k) = VY;(k,, k^, k,,) + UYi,i(k„ k„ k,), (9b) 

24 where i=l,2,...,N-l, and where ki = 1,2,. ..,M, and where VX,=X(q-N+l), 
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1 VY,=Y(q-N+l) are quadrature components of the first received symbol in the sequence 

2 of N processing symbols. 
3 

4 A final N-symbol decision regarding the q-th transmitted symbol corresponds to 

5 the maximum of a calculated vector at the i=(N-l)-th step of the recurrent procedure, i.e. 
6 

7 Max {[VX,,(k„ k^, . . K^Of + [VY^Ck,, k^, . . k^ J]'}, (9c) 

8 

9 where the maximum is determined through all "k". 
10 

11 It should be noted that decision algorithm (9c), which is similar to algorithm (6c) 



1 2 for 3-symbol processing, may be used for both a symbol-by-symbol decision and a joint 

13 decision on (N-1) symbols simultaneously. In the case of a symbol-by-symbol decision, 

14 maximization in (9c) is provided through indexes k^.j, corresponding to phase differences 

1 5 between the (N-l)-th and N-th symbols. In the case of the joint symbol decisions, 

1 6 maximization in (8c) is provided through the whole set of indexes k,, kj, . . . , k^.i, 

17 corresponding to phase differences between all N symbols participating in N-symbol 

1 8 processing procedure. 
19 



20 To reiterate, according to the invention, an N-symbol noncoherent algorithm for 

21 DMPSK signals can be described as follows: 
22 

23 a) N consecutive symbols with indexes q, (q-1), . . ., (q-N+1) participate in making a 

24 decision about the current q-th symbol, and the corresponding quadrature 

25 components of the received signals X(q),Y(q); X(q-l),Y(q-l); . . .; 

26 X(q-N+l),Y(q-N+l) are saved; 
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1 b) A set of quadrature components X(q),Y(q); X(q-l),Y(q-l); . . .; X(q-N+1), 

2 Y(q-N+ 1 ) are subjected to a (N- 1 )-step transformation; 

3 c) The i-th step, where i= 1 ,2, . . . ,N- 1 , includes 

4 ( 1 ) Calculation of multi-dimensional transforms UXi+,(k, , kj, . . . , k^) and 

5 UYi+i(k|, kj, . . ki) of the current received components X(q-N+i+l) and 

6 Y(q-N+i-i-l) according to equations (8a) and (8b); and 

7 (2) Calculation of mult-dimensional components VXi+,(k,, k2,. . k^) and 

8 VYi^.,(ki, kj,. . kj) according to recurrent formulas (9a) and (9b); 

9 d) Final decision relative to the current q-th symbol is the phase difference from 

10 the set of differences A(kN.i), where kN.,=l, 2,. . M, corresponding to the 

1 1 maximum of the length of vectors [VXN(k,, kj, . . k^.i), VYN(ki, kj, . . k^.i)], 

1 2 according to algorithm (9c); 
13 

14 Turning now to Fig. 6, a schematic diagram is provided of the N-symbol 



1 5 noncoherent demodulator of DMPSK signals which implements equations (9a) and (9b). 

1 6 The apparatus of Fig. 4 includes (N-1) Intersymbol Processors IPj . . .IPn-i; (N-2) pairs of 

1 7 memory blocks VX2,VY2 . . • VXn.i,VYn.i, a Decision block, and two memories for signal 

18 quadrature components X and Y, which are depicted as Shift Registers SR, and SRy. 
19 

20 The apparatus of Fig. 6 operates as follows. Quadrature components of the 

21 current received symbol X(q) and Y(q) are directly fed to the IP^.!. In addition, 

22 component X(q) is inserted into the first (right-hand) cell of SR^, shifting all register 

23 elements by one cell from right to left, and component Y(q) is inserted into the first 

24 (right-hand) cell of SRy, shifting all register elements by one cell from right to left. As 

25 seen in Fig. 6, storage locations of both memory shift registers SR^^ and SRy are connected 

26 to corresponding intersymbol processors IP. Each IP, in turn, receives signals from 

- 24 - 



PCTEL-025 



1 registers SR^ and SRy as well as results of computations of the previous IP stored in 

2 corresponding memories VXj and VYj. All intersymbol processors provide computations 

3 according to equations (8) and (9a) and (9b}, but the number of operations carried out by 

4 each IP is different. For example, IP, calculates M complex numbers, and IP^., calculates 

5 M*^^ complex numbers per symbol. Thus, for 4-symbol processing of DQPSK, the first 

6 IP will calculate four numbers, while the last IP will calculate sixty-four numbers. As 

7 shown above, in the proposed algorithms all computations are simple summing 

8 operations of two numbers (see Tables 1-5). A final decision is provided in the Decision 

9 block according to algorithm (9c). 
10 

11 It should be appreciated by those skilled in the art that Fig. 6 illustrates only the 

1 2 functions of the proposed multisymbol noncoherent demodulator and not necessarily its 

1 3 most desirable implementation. Indeed, there are numerous ways to implement the 

14 proposed recurrent algorithm of equations (9a) and (9b) and the sequence of operations 

1 5 provided by apparatus in Fig. 4. For example, all operations provided by the individual 

1 6 intersymbol processors of Fig. 6 can be combined in a single high-speed intersymbol 

1 7 processor. Likewise, all memory and shift register storage can be provided in a single 

1 8 block of memory. 
19 

20 Fig. 7 provides a high level block diagram of an alternative implementation of the 

21 multisymbol demodulator of the invention. The apparatus of Fig. 7 includes a single 

22 intersymbol processor IP (which may be implemented in a DSP, discrete components, or 

23 in software of a general-purpose processor), X-memory and Y-memory blocks (which 

24 may be part of a single memory such as a RAM, or separate memory elements), and a 

25 Decision block. The X-memory stores components X(q-N+2), X(q-N+3), . . ., X(q-l) of 

26 the received signal as well as intermediate results of intersymbol processing VXj 
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1 (j=2,3, . . . ,N- 1 ). Similarly, the Y-memory stores components Y(q-N+2), Y(q-N+3), . . . , 

2 Y(q-l) of the received signal as well as intermediate results of intersymbol processing VYj 

3 (j=2,3,...,N-l). The Intersymbol Processor (IP) calculates all intermediate components 

4 VXj and VYj as well as final values VX^ and VY^ according to equations (8) and (9a) and 

5 (9b). 
6 

7 The DBPSK and DQPSK modulation techniques described above were simulated 

8 by means of stochastic simulation programs which were progranuned to compare perfect 

9 coherent processing with noncoherent processing having different processing intervals 

10 according to the invention. Simulation results for DBPSK and DQPSK modulation 

1 1 techniques are shown in Fig. 8 and Fig. 9. More particularly. Fig. 8 shows the bit error 

12 rate (BER) as a function of signal-to-noise ratio (SNR) in dB for DBPSK modulation. 

1 3 The lowest (solid) curve corresponds to perfect coherent processing of the DBPSK signal, 

14 while the remaining curves (from top to bottom) correspond to noncoherent processing 

15 according to the invention with 2, 3, 4, and 5-symbol intervals. 
16 

17 At least two conclusions can be deduced from the simulation results. First, in the 

1 8 case of DBPSK there does not appear to be any reason to increase the processing interval 

19 beyond five symbols, because with 5-symbol interval the system already approximates 

20 minimum possible BER, which is approximately twice the BER for the perfect coherent 

21 processing of BPSK. This conclusion is illustrated in Table 8, which shows BER vs. SNR 

22 for BPSK with perfect coherent processing, DBPSK with perfect coherent processing, and 

23 DBPSK with 5-symbol noncoherent processing. 
24 

25 
26 
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1 Table 8 BER vs. SNR for BPSK and DBPSK 
2 



SNRdB 


BER BPSK 

(Perfect Coherent) 


BER DBPSK 
(5-symbol noncoherent) 


BER DBPSK 

(Perfect Coherent) 


7.5 


4.0 E-4 


1.0 E-3 


8.0 E-4 


8.0 


2.0 E-4 


4.5 E-4 


4.0 E-4 


8.5 


8.5 E-5 


1.8 E-4 


1.7 E-4 



3 

4 A second conclusion is that in the case of DBPSK the most considerable 

5 improvement takes place at the transition from 2-symbol processing to 3-symbol 

6 processing; it being evident from Fig. 8 that the dashed line of 2-symbol DBPSK is 

7 further away from the dashed line of 3-symbol DBPSK than the dashed line of 3-symbol 

8 DBPSK is from the dotted line of 4-symbol DBPSK. For example, at a BER of 4x10"^, 

9 going from 2-symbol DBPSK to 3-symbol DBPSK yields a gain of slighdy over .3dB, 

10 whereas increasing from 3-symbol DBPSK even to 5-symbol DBPSK will yield a gain of 

1 1 only about another . IdB. So, from a practical point of view, it may be desirable to use 3- 

12 symbol pilotless noncoherent processing for best performance gain with limited signal 

13 processing costs. 
14 

15 Fig. 9 shows the bit error rate (BER) as a function of signal-to-noise ratio (SNR) 

16 in dB for DQPSK modulation. The lowest solid curve corresponds to perfect coherent 

17 processing of the DQPSK signal, with the remaining curves (from top to bottom) 

18 corresponding to noncoherent DQPSK processing with 2, 3, 4, and 5-symbol intervals. 

19 As was the case for DBPSK, the most considerable improvement in DQPSK takes place 

20 when going from 2-symbol processing to 3-symbol processing. Further improvement 

21 decreases as the number of processing symbols increases. However, practically, 4- 
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1 symbol processing is a good option for DQPSK modulation, because the gain is still 

2 rather significant between 3-symbol and 4-symboi processing. For example, at a BER of 

3 7x10'^, going from 2-symbol processing (the top dashed line) to 3-symbol processing (the 

4 next dashed line) yields a gain of almost .5dB. Going from 3-symbol processing (the 

5 second dashed line) to 4-symbol processing (the dotted line) yields a similar, albeit 

6 slightly smaller gain of over .4dB. However, going from 4-symbol processing to 5- 

7 symbol processing yields only approximately .IdB. Thus, 4-symbol DQPSK would 

8 appear to give the best performance gain relative to the additional complexity in 

9 processing. 
10 

1 1 There have been described and illustrated herein several embodiments of 

12 methods, apparatus, and systems for implementing non-coherent processing of pilotless 

13 telecommunication signals. While particular embodiments of the invention have been 

14 described, it is not intended that the invention be limited thereto, as it is intended that the 

1 5 invention be as broad in scope as the art will allow and that the specification be read 

16 likewise. Thus, while particular hardware and logic elements have been disclosed for 

17 implementing the invention, it will be appreciated that the invention could be 

18 implemented using different hardware and logic elements as well. In addition, the 

19 invention can be implemented with software and firmware as well. It will therefore be 

20 appreciated by those skilled in the art that yet other modifications could be made to the 

21 provided invention without deviating from its spirit and scope as claimed. 
22 
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